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Abstract. We consider the implications of the condensation of a general local BRST invari- 
ant dimension two operator built out of the localizing ghost fields of the Gribov-Zwanziger 
Lagrangian which is a localized Lagrangian incorporating the Gribov problem in the Landau 
gauge. For different colour tensor projections of the general operator, the properties of a frozen 
gluon propagator and unenhanced Faddeev-Popov ghost propagator, which are observed in lat- 
tice computations, can be reproduced. The alternative possibilities are distinguished by the 
infrared structure of the propagators of the spin-1 fields, other than those of the gluon and 
Faddeev-Popov ghost, for which there is no numerical simulation data to compare with yet. 



1 



1 Introduction. 



Yang- Mills theories play an important role in our understanding of the fundamental particles of 
nature. For instance, the strong force is believed to be reliably described by an SU (3) non-abelian 
gauge theory whose fundamental particles are quarks and gluons. In nature these particles 
cannot be isolated as separate entities. Instead at high energies they are asymptotically free, 
[HE], but at low energies they are present only within bound states of hadrons. Understanding 
the specific dynamics of the mechanism which confines quarks and gluons and prevents them 
from being observed in nature is currently a major topic of interest. Ordinarily in a quantum 
field theory the behaviour of a particle's propagator plays a key role in its interpretation as an 
observed quanta. For instance, in Quantum Electrodynamics (QED) the electron propagator has 
a simple pole at real positive value of its squared momentum which corresponds to its physical 
mass when all quantum corrections have been computed to all orders in perturbation theory. 
The form of this fundamental propagator is therefore associated with a real observed particle. 
By contrast in Quantum Chromodynamics (QCD), which is the non-abelian generalization for 
the strong force, the propagators of the quark and gluon fields derived from the canonical 
Lagrangian emerge as being of the same fundamental form but in the case of the gluon it is 
massless. Clearly this cannot be correct since that would imply that the gluon should exist as 
a free massless observed particle. However, this analysis derives from a gauge fixed Lagrangian 
with a non-abelian symmetry. It transpires that contrary to what occurs in an abelian theory 
there are problems in fixing the gauge uniquely globally in Yang-Mills theories. This was first 
pointed out in detail by Gribov in his seminal work, [3]. In essence one can only fix a gauge 
uniquely locally in a non-abelian gauge theory but not globally. It turns out that one can 
construct different gauge configurations satisfying the same gauge condition in the non-abelian 
case. Therefore, there is an overcounting in the definition of the path integral and the region of 
integration in configuration space has to be restricted to a specific subspace. This is known as 
the first Gribov region and it contains the origin, [3]. Whilst this overcomes a significant amount 
of the copy issue the gauge is still not uniquely fixed. Only inside the fundamental modular 
region which is within this horizon is there a unique globally fixed gauge. 

In analysing this problem for the Landau gauge Gribov, [3j, managed to determine several 
interesting features. First, the path integral could be modified to incorporate the restriction to 
the Gribov region. The boundary is defined by the no-pole condition where the Faddeev-Popov 
operator vanishes. This means that the inverse is finite within the horizon of the Gribov region. 
Using a semi-classical analysis Gribov demonstrated that the path integral cut-off modified the 
gluon propagator. In particular the fundamental behaviour was replaced by a propagator which 
vanished at zero momentum, known as gluon suppression, but the simple massless pole property 
accepted at high energy was retained. Underlying this was a new mass parameter called the 
Gribov mass and denoted by 7, [3]. It is not an independent quantity but is a function of the 
coupling constant, g, and defined via the horizon condition cutting off the path integral leading to 
a gap equation. Only when the Gribov mass actually satisfies the gap equation can the theory 
be regarded as a gauge theory. This gap equation was responsible for another novel feature 
which is that the Faddeev-Popov ghost propagator behaved as a dipole in the infrared which 
is known as ghost enhancement. This latter property was subsequently reformulated in terms 
of the Kugo-Ojima confinement criterion, [HE], which was established by a BRST analysis of 
the Landau gauge. However, one of the main consequences is that the analysis showed that the 
gluon propagator ceased to be fundamental when the copy problem was taken into consideration 
and the hope was that the Gribov problem was fundamental to the problem of confinement, [3]. 

One of the main difficulties in extending the semi-classical analysis of [3] was that the no pole 
condition introduces a non-local operator into the Lagrangian. This is an obstacle to performing 
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canonical field theory calculations which only proceed when a Lagrangian is at least local and 
renormalizable. However, in a series of articles, [61 [TJ [U [91 [TUl [UJ 121 [131 CG3 HE], Zwanziger 
overcame this by managing to not only rewrite the Gribov Lagrangian in a local form but also in 
a way that it was renormalizable, [13 11171 [T8]. The localization proceeded by the introduction of 
a set of additional ghost fields, {4>^, <^f, u^ b , u)" t 6 }. The first pair are bosonic whilst the second 
set have Fermi statistics and are necessary to ensure the theory remains asymptotically free 
at high energy, for instance. The renormalization group functions are not affected by these 
additional fields whose effect is only apparent in the infrared region, [13 1 117 1 ITS]- Moreover, no 
extra renormalization constants are required since the anomalous dimensions of the localizing 
ghosts are the same as that of the Faddeev-Popov ghost and 7 is rendered finite merely by a 
simple combination of the gluon and ghost renormalization constants. With these additions to 
the Lagrangian the original suppressed gluon propagator is retained. However, one can study 
the loop corrections using the localized Gribov version which is known as the Gribov- Zwanziger 
Lagrangian. For instance, the one loop gap equation can be computed, [8], [T3] , and reproduces 
that of Gribov, [3]. This was extended to two loops in the MS scheme in [19] with the one loop 
propagator corrections also deduced in [20l [21] . As noted in [HI [12] by Zwanziger the gluon 
remains suppressed when quantum corrections are included. Also the Kugo-Ojima criterion was 
shown to be satisfied at two loops, [T9l [20] . Subsequently this criterion was revisited but in the 
context of the Gribov-Zwanziger Lagrangian, [22J, where the additional localizing fields have to 
be included in repeating the original BRST analysis of [U[5]. It was shown that the Kugo-Ojima 
condition was valid in the Gribov-Zwanziger Lagrangian. 

Having summarized the theoretical analysis of accommodating the gauge copy problem in 
the Landau gauge in QCD the next issue in this area relates to lattice data collected over 
the last few years on the gluon and Faddeev-Popov ghost propagators. Given the advances in 
computing power and algorithms to carry out the gauge fixing numerically on larger lattices, 
[231 [Ml 1251 1261 [271 1281 1291 [301 131] , there is now a reasonable amount of data on both propagators 
at low momenta. More specifically, the zero momentum behaviour appears to show clearly 
that contrary to the Gribov scenario the gluon propagator freezes to a non-zero value and the 
Faddeev-Popov ghost does not enhance. Instead its infrared form does not significantly deviate 
from that of its ultraviolet fundamental massless form. Indeed there is also similar evidence from 
Schwinger-Dyson computations. See, for example, [311 1321 [33] . Thus there are two distinct cases. 
One is the gluon suppressed but ghost enhanced propagators known as the scaling or conformal 
solution. The other is the non-zero frozen gluon propagator with the unenhanced ghost. Such a 
gluon propagator was originally derived in [33] in a Schwinger-Dyson analysis and referred to as 
the massive solution. More recently, [34] . the term decoupling solution has been used. As both 
labels are recognised we will use them synonymously in the discussion. Therefore, there has 
been a debate as to how the numerical data fit in with the Gribov scenario, [3], and whether the 
absence of the Kugo-Ojima confinement was significant. One point of view was that satisfying 
the criterion or not was an additional condition on the gauge fixing procedure, [35], effective in 
the infrared. However, given the properties and elegance of the Gribov-Zwanziger Lagrangian 
it seemed appropriate to try and extract a gluon propagator with non-zero freezing and a non- 
enhanced Faddeev-Popov ghost in that approach. This was achieved in the series of articles, 
[361 [37l [38] . In essence the Lagrangian was refined to include a local BRST invariant dimension 
two operator built from the localizing ghost fields. Ordinarily such a term would just introduce 
masses for these fields only. However, the complicated way the localization arises means that 
in constructing the propagators from the quadratic sector the extra mass term affects the gluon 
propagator too. Specifically the extra mass parameter appears in the gluon propagator in such 
as way that there is non-zero freezing at zero momentum. Equally the extra mass excludes the 
enhancement of the Faddeev-Popov ghost and therefore the numerical data and Schwinger-Dyson 
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analyses can be modelled by what was termed a refinement of the original Gribov-Zwanziger 
Lagrangian. The justification for the inclusion of this additional mass operator is to consider 
the condensation of the operator via a dynamical mechanism. The analysis for this is achieved 
by the local composite operator (LCO) formalism, |39|, I40j . which was initially used in QCD to 
study the condensation of a gluon mass operator, [4T| [42]. 

Whilst the computation of [36], [37l [38] clearly accommodates the current picture of the gluon 
and ghost propagators it transpires that the most general dimension two BRST operator was 
not considered. Given that the localizing fields carry more than one colour adjoint index there 
are six different colour projections. Therefore, it is the purpose of this article to perform a 
comprehensive analysis of the various different colour channels. It will turn out that there are 
other operators whose condensation will lead to a non-zero freezing of the gluon propagator 
and an unenhanced Faddeev-Popov ghost propagator. Therefore, the refined Gribov-Zwanziger 
Lagrangian of [361 [371 [38] is not the unique model of the lattice data which is why we will refer to 
the explanations here as an alternative. Moreover, in our results we will show that the conformal 
or scaling behaviour of the original or pure Gribov-Zwanziger analysis can persist even with the 
extra mass operator. Although the alternative we focus on here, which will be referred to as the 
7Z channel for reasons which will be clear later, is equivalent to the earlier explanation, which will 
be called the Q channel, the full infrared behaviour of both are not equivalent. For instance, the 
behaviour of the Bose ghost propagator is different in each channel. The study of this particular 
propagator has been of interest recently in the pure Gribov-Zwanziger Lagrangian, |43[ I21j. It 
has been shown that there is enhancement in certain colour channels. This was observed at the 
one loop level, [21] . but proved to be an all orders feature in [43] , Briefly, the enhanced Bose 
fields are the Goldstone bosons associated with the spontaneous breaking of the BRST symmetry 
in a theory where the fields are constrained by the horizon condition, [33]. In principle one could 
resolve which of the Q or 1Z cases was correct when the BRST invariant operator is included if 
there was numerical data for this propagator. However, these localizing ghost fields only exist 
within the Gribov-Zwanziger Lagrangian itself and not within the lattice formulation of QCD. 
Therefore, it is not clear whether one could construct the equivalent correlation function in 
the lattice language. Whilst we will carry out the one loop analysis for the 1Z channel we will 
also give a general indication of the properties of the additional colour channel possibilities in 
order to have as a complete picture as possible. This would become important if additional 
numerical data became available which rules out either of the Q or 1Z channel explanations. 
Interestingly in certain scenarios the colour tensor structure which was observed in |43[ 12 lj , but 
for the enhanced Bose localizing propagator at one loop, actually arises in the dominant term 
in the zero momentum limit in the original propagator prior to any loop analysis. 

The article is organised as follows. The key properties of the pure Gribov-Zwanziger La- 
grangian are reviewed in the next section. The most general local BRST invariant dimension 
two operator built from the localizing ghosts is introduced in section three and the propagators 
are constructed for each of the six potential single additional mass parameters. The justification 
for the inclusion of such additional masses is provided in the subsequent section where the focus 
is on the 7Z channel. There the one loop effective potential for the operator is constructed as 
well as the one loop MS corrections to all the propagators. These are required to study the 
effect of the gap equation satisfied by 7 has on the propagator corrections in order to see, for 
example, whether there is any Bose ghost enhancement. Our conclusions are given in section 
five. There are two appendices. The first gives the full propagators for the inclusion of all six 
colour channel projections of the operator for the case of SU(3) whilst the second records the 
propagators for the W and S channels separately for SU(N C ). Whilst the expressions recorded 
in both appendices for the propagators are cumbersome they are significantly compact compared 
with those for the general colour group case. 
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2 Background. 



We begin by briefly recalling the essential features of the Gribov-Zwanziger Lagrangian and the 
construction of the propagators. Throughout the article we work in the Landau gauge. The 
basic Lagrangian is, [3], 

L Grib = L QCD Crf 1 _ dN^ 

2 d v D u * 2g 2 y 1 

where the original gauge fixed Lagrangian is 

LQ CD = - ^G a ^ u G a ^ - _L(^AJ) 2 - c a d f "D fl c a + i^p^ 1 (2.2) 

and a is the covariant gauge fixing parameter. Although we have included a we will focus 
purely on the Landau gauge, corresponding to a = 0, but it is retained here since it is crucial 
in obtaining the propagators. The parameter 7 is the Gribov mass which is not an independent 
quantity but is the mass scale which is implicit in the horizon condition, [3], 

where is the dimension of the adjoint representation. The presence of the horizon condi- 
tion alters the infrared behaviour of the gluon and Faddeev-Popov propagators in the original 
Gribov scenario, [3]- As has already been well documented the non-local term renders practical 
computations impossible. However, Zwanziger's series of articles on how one can localize the 
non-locality produces a renormalizable local Lagrangian, [3 EJ [HJ H21 131 HH SI] , which is 

L GZ = L QCD + l p *b»d v (D vPft ) ab + ~p abti d v (D v ^) ab ~ l -e b ^d u {D vPll ) ab 

- ±=gf^ d u^e {DvC) b e c, _ i7 2 f abc A a^bc _ <£^_ {2A) 

There are additional fields to the gluon, A®, Faddeev-Popov ghost, c a , and massless quark, ip l1 , 
which are two Bose ghosts, and p" 6 , and the Grassmann ghosts lu^ and u'jf. The indices 
range over the values 1 < a < Na, 1 < I < Np and 1 < i < Nf where Np is the dimension of 
the fundamental representation and Nf is the number of massless quarks. In earlier localized 
versions of the Gribov Lagrangian the complexified Bose ghosts 4> ab and (p^ were used. Given 
recent developments, [431 44J, in terms of the propagator structure of the real versions of the 
Bose ghosts we will use this formulation but note that the relation between the two versions are 
trivial and given by 
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^ (pf + <) , 4>f = ^ ~ <) • (2-5) 



However, the renormalizability proofs of the Lagrangian, [131 [T71 118j . were performed for the 
complex Bose ghost field and used the BRST symmetry of the localized Lagrangian. Since we will 
be considering a BRST dimension two operator we note the BRST symmetry transformations 
for flZH) given (|23]> are 



SA% = - (D„c) a , 5c a = ^f abc c b c c , 5c a = b a , 5b a = 

5<t>f = < , 5ujf = , 5^ b = , 5ojf = 4>f (2-6) 
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where b a is the Nakanishi-Lautrup auxiliary field. Therefore, it is trivial to observe that the 
colour non-singlet operator 

Qabcd 



u ab u cd 



4> ah 4) cd 

is BRST invariant. Rewriting in terms of the real Bose ghosts the operator is equivalent to 

■\abcd 



2.7) 



0° 



p ab p cd + ^ab p cd _ if) ab^cd + ^ab^cd 



uj ab uj cd 



2.8) 



In |36t [37] [38] it was the condensation of the operator § ac (j bd Q abcd which was investigated using 
the local composite operator formalism in order to see what effect it had on the infrared structure 
of the propagators. 

In this localized Lagrangian the horizon condition, (|2.3p . is replaced by the vacuum expec- 
tation value of local fields, [EJ IH], 



idN A j 2 

which determines 7 to two loops in the MS scheme, [19], by the solution of 
[5 3 (C A ^W 

[ 35 



(2.9) 



1 



r 2 



3893 22275 29 . > 65 , / C A j 4 

-s 2 H C(2) In / 

1536 4096 128 w 48 \ ^ 

411 f- 1317^ 2 \ 



+ 



128 



In 



C A 1 A 



1024 



4096 



, W |^-l-ce> + ^ 74 ' 



I In 



+ 0(a 3 ) 



(2.10) 



where Q{z) is the Riemann zeta function, S2 = (2\/3/9)Cl2(27r/3) and Cl2(x) is the Clausen 
function and a = g 2 /(16ir 2 ). However, it was demonstrated recently that the horizon condition 
could be replaced by a vacuum expectation value of an operator involving only the Bose ghost, 
albeit one involving an infinite number of terms, [25]. Specifically, 



jabcd / Qu^abfi 



dC A N Al A 



ch 4 

+ 0(g 3 ) 



r 



(2.11) 



where f£ bcd = jo-bpjcdp^ 'pj^^ reproduces (|2.10p to two loops and is constructed by application 
of the equation of motion 

\f ahc {d v D v ^) hc (2.12) 



4" 



Cat 2 



to ()2.9p via the intermediate vacuum expectation value 

f ab Pf cd P ^ab, {x) (PD^cd^ 



dC A N A1 A 
9 2 



(2.13) 
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where d is the spacetime dimension. 

One of the reasons for reviewing the various formulations of the gap equation is to motivate 
our alternative mechanism for modelling the massive or decoupling solution observed on the 
lattice. First, we recall that the Landau gauge propagators for the spin-1 fields of (12. 4p are 

5 ° V *Wp> 





-p)) = 




-p)> = 


(A«(p)p b v c (-p)) = 




-p)> = 




-p)> = 


(pf(p)p C u(- 


-p)> = 



[b 2 ) 2 + c A7 4 

z/ a6c 7 2 



[Cp 2 ) 2 + C A7 T^ 




p 2 P 2 [(p 2 ) 2 + Ca7 4 ] 







where 



= - ^-17^ (2-14) 

P^(p) = ^ ~ ^ • V(P) = ^ (2-15) 

are the usual transverse and longitudinal projectors. Therefore to reproduce the leading order 
contribution to the gap equation, and hence obtain Gribov's original expression, one simply 
integrates the mixed A a ^-£^ propagator of (|2.14p . Equally the alternative formulation of the 
gap equation involving only the Bose ghost can be justified at leading order by noting that the 
second term of the propagator needs the massless pole to be absent. This is achieved by 
the inclusion of the wave operator, d a d a , at leading order, [3], and eventually correctly with the 
Faddeev-Popov operator in keeping with the ethos of Zwanziger's extension of Gribov's semi- 
classical horizon condition, [7j. Given this, one can obtain an idea of which of the other operators 
may condense and be dominant or relevant in the infrared at leading order by considering other 
propagators of (|2.14p . Clearly integrating the gluon propagator suggests the condensation of the 
gauge variant operator ^A^A a ^ which has been widely treated. However, if one considers the 
BRST invariant operator, (12. 8ft . then integrating the relevant Bose and Grassmann propagators 
of (|2.14p we find the colour dependence 

(O abcd ) oc f abe f cde . (2.16) 

Clearly the other colour channel of the relevant propagators, 5 ac 5 bd , gives zero upon integration 
in dimensional regularization. Whilst we are not saying that there could be no condensation 
in this colour channel, or any others not present here, we are suggesting that the more natural 
colour channel to condense is that of (|2.16p . It is natural in the sense that it does not differ 
too much from the original localized Lagrangian structure. Moreover, it would appear to be the 
dominant expectation value at leading order. Given this simple observation we will examine the 
consequences for this taking this colour channel choice in (|2.8p and discuss alternative cases as 
well in the next section. 



3 Propagators. 

In this section we consider the construction of the spin-1 propagators for (|2.4p when there are 
extra mass terms eminating from a BRST invariant mass. As Q abcd is BRST invariant then in 
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order to have a colour scalar it needs to be contracted with a rank 4 colour tensor. To be as 
flexible as possible at the outset we take the general operator 



O 



^ 5 ac 5 bd + ^ f ace f bde + f^ f abe f cde + ^abcd + f^gobgcd + ^gadgbc 

Ca Na 



Qabcd 



(3.1) 

where /if are the various mass parameters. The labelling is chosen in order to track the location 
of the various colour structures within the final propagators. The inclusion of Ca and Na in 
several terms is in order to have a degree of uniformity in the various different final propagator 
forms and to make relative comparisons easy to follow as will be evident later. In addition to 
be complete we include an additional dimension two operator which is 

A 2 = ^A a ^A a » . (3.2) 

In a general linear covariant gauge for (|2.2|) the operator ^A°^A ail — ac a c a is BRST invariant 
and was introduced as a potential gluon mass operator by Curci and Ferrari in [36]. So in 
the Landau gauge (13.2j) represents a natural additional operator in the context of studying 
the BRST operator for the localizing ghost sector. So it seems appropriate to include it in 
our current construction to ascertain its effect. The subscripts, X, Q, W, TZ, 5, V and T, 
derive from notation used in earlier papers, [201 [21]. There the colour structure of the one loop 
corrections to the 2-point functions was examined and the inversion of the appropriate matrix 
of 2-point functions was performed to obtain the one loop propagator corrections. This involved 
understanding the multiplication of the colour tensors associated with each amplitude. As the 
colour structure for the BRST invariant mass term is now of the most general form, we will use 
the same approach as [20, 21J to construct the propagators and recall the essentials for this. 



First, for the spin-1 sector we focus on the fields which mix through the term involving 7 2 in 
(|2.4jl . Therefore, we take our basis of fields to be {A^, p ab } in that order. Although there is 
no mixing with p ab we have included it here since there is potentially an effect in the longitudinal 
sector. Next we define the matrix of quadratic terms in the momentum space version of the 
Lagrangian as, [20] . 

x5 ac U pcd o 

A {ab\cd} = | U fcab Qabcd q | (33) 



Q p 



abed 



for the transverse sector and 



(p^Lgac jj^L jacd jacd \ 

U L f cab QLabcd q ^3^ 

yL jcab q qL abed J 

for the longitudinal sector where we will use the superscript L throughout to differentiate it 
from the transverse sector. The respective Lorentz projectors are passive and omitted, |20|, I21j. 
The colour decompositions are taken to be 

Q abcd = Q^ac § bd + yy^acejbde + p^abe jede + globed + p^abged + p^dgbc 
Qabcd = Q p5 ac 5 bd + ^jacejbde + ^jabejede + globed + p p§ ob § cd + J^gadtfc 

(3.5) 

which is the origin of our earlier notation and 

■\Labcd s^iLxacxbd , -\a;L face fbde , tjL fabe fede 1 eLjabcd , TjLxabsrcd , TLzadxbc 



QLabca QLgacgbd _|_ yyL jace jbde _|_ pj^ jabe jcae _|_ gL jabed _|_ spLgabgcd _|_ pL^adgb 

5 b 

(3.6) 



QLabcd QLgdcgbd _|_ yyL jace jbde _|_ pL jabe jede _|_ gL^jdbcd _|_ spL gab gcd _|_ pLgadgbc 



s 



with 



jabcd 

a A 



±Tr (T%T%T%T d 2) 



(3.7) 



being the rank 4 totally symmetric tensor in the adjoint representation, [47J. The propagators are 
then obtained by inverting the quadratic part of the momentum space Lagrangian. However, this 
is more involved than usual due to the colour structure and is formally given by the corresponding 
colour structures for the transverse sector, 



jj{cd\pq} 



A6 cp gfcpq 
jpcd jycdpq 






jycdpq 



(3.8) 



and 



■qL {cd\pq} 



( A L 5 cp B L f cpq C L f cpq 

jpcd jyL cdpq gL cdpq 
^ qL jpcd gLcdpq 



(3.9) 



cL cdpq 



for the longitudinal part of the propagators. The two sectors of the propagators can be split 
since we employ the projectors P pu (p) and L pu (p) which satisfy 



V = P Hv(p) + L pv{P) > P ^ip) L vaip) 

We use a similar decomposition for the colour tensors with 



. 



(3.10) 



T, cd P1 
jycdpq 



V^5 Cp 5 dq + J^efqe + ^jcde fpqe + L ^cdpq + ^cdgpq + tf^cqgdp 

V p 8 cp 5 dq + j p f cpe f dqe + )C p f cde f pqe + £ p d% pq + M p 8 cd 5 pq + M p 5 cq 8 dp 



and 



sp.L cdpq 

u i 

cL cdpq 
cL cdpq 



(3.11) 



V L 6 cp 5 dq + jLjcpejdqe + jcde jpqe + jT^cdpq + M L § cd § pq + ^f^cq^p 



cL 



+ T^f Cpe f dqe + gLfcdejpqe + ^cdpq + yLgcdgpq + fgL 
F Ljcpe f dqe + gL f cde jpqe + % L ^cdpq + yLgcdgpq + 



The inversions in the two sectors proceed via 



^{ab\cd}j^{cd\pq} 




and 



j^L {ab\cd}j^L {cd\pq} 



5 cp 
5 cp 5 dq 
5 cp 8 dq 



(3.12) 



(3.13) 



(3.14) 



where the matrix on the right hand side is the unit matrix in this colour space basis. In order 
to handle the products of the colour tensors in this matrix multiplication we recall, [21], that 



,abpq,cdpq 

a A a A 



ai§ ab 6 cd + a2 f S ac S bd + 6 ad 6 bc 



f 



a pa 



bl6 ab § cd + b l § ac § oa + § aa § 



iac rbd 



cad <;bc 



| ^ / jace jbde _|_ jade jba 



+ b 4 d a A bcd 



a 3 (f ace f bde + f ade f bce ^j + a 4 d a A bcd 

(3.15) 
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where 



ax 



hAQC 2 A N A {N A - 3)d a A bcd d A pq d a A bpq + U4{2N A + 19) (df cd d a A bcd y 



150C A N A (3N A + ll)d a A bcd d a A bcd + Q2bC\N'i 



MN A (N A 



3)[l2(N A + 2)d e A f9h d e A f9h 



25C A N A ] 



0-2 



144(117V A - 8) (d a A bcd d a A bcd ) - 1080C 2 A N A (N A - 3)d a A bcd d cdpq d 



+ 625C A Nj - 3000C A N A d A bcd df 

1 



•abed jedpq pbpq 



W8N A (N A - 3)[12(N A + 2)d A f9h d A fgh - 2hC\N A ] 



'abed jabed 



«3 



a i 



[12(N A + 2)d a A bcd d 



25C\N A ] 



5AC A N A (N A - 3) 
[216(N A + 2)d A bcd d A pq d A m - 125C A N A - 360C 2 A d a A bcd d a A bcd } 
18[12(iV A + 2)d\ S9h d e l qh - 25C\N A ] 



and 



[5C\N A - 12d a A bcd d a A bcdi 
9C A N A (N A - 3) 



2b 2 



[6(Na 



l)d a A bcd d a A bcd - 5C A N A ] 



9C 2 A N A (N A - 3) 



(3.16) 



- 9a 

3 

(3.17) 



It is clear that the set of linear algebraic equations resulting in multiplying out the matrices 
of colour amplitudes will be very involved. For reference in the case of the absence of any 
conventional mass terms these are formally given in [21\ 145] . Indeed retaining all possible 
masses /x 2 will be very complicated and we have provided the propagators for this situation in 
Appendix A for the specific case of SU(3). Instead it seems more instructive to consider the 
effect one particular mass term has on the propagators in turn. The motivation for this is to see 
which masses can produce propagator behaviour akin to that observed in lattice simulations. 
We note that in solving the set for the longitudinal sector it is important one follows a specific 
algorithm. This is because we are working in the Landau gauge but in order to carry out the 
inversion correctly we must retain a non-zero a at the outset. The Landau gauge propagators 
are deduced at the end by setting a = which will produce a transverse gluon propagator in 
all cases. As an aid we note that in the absence of (I3.ip and (|3.2|) the non-zero propagators are 



(A»(p)4(-p)> 

<^»e(-p)> 
(A^p)p b u c (-p)) 

(C b (p)pl d (-P)) 



abjl 



6 ab p 



ab^2 



if abc J 2 



a5 p 



[(p 2 ) 2 + C7 A7 4 ] 



Pauip) + 



[(p^ + aC A ^] 
iaf abc -f 2 



L^(p) 



[(p^ + aC Al ^ {p) 







acxbd 



jabe jcde^yA 



a jabe jcde^4 



P 2 [(p 2 ) 2 + C A ^y^ yFJ " p 2 [(p 2 ) 2 + aC A 7 41 



(pt(p)pl a (-p)) 





(u ab (p)ul d {-p)) 



aczbd 



5 ac 5 



(3.18) 



which are clearly non-singular in the limit to the Landau gauge. 



10 



Given these considerations we record the propagators for each of the masses pq being non- 
zero in succession. We append a subscript i to the propagators themselves to keep a track of 
each channel. First, before looking at the six localizing ghost possibilities, including the gluon 
mass we find 

5 ab V 2 

(A«(P)AK-P))x = ~ [{p2)2 + ^ p2 + CA ^ (p) 

Wp)#(-p)>* = WWT^ P ^ 

{Al(p)p b v c (-p)) x = 

xacxbd fabe fcde^.4 

(Q b (p)t?(-P))x = — + w ( J, 2 2, r 4T P ^(P) 

p p a p z [(p z y + h z x p + Cai \ 

xaczbd 

(P; b (p)pl d (-P))x = «(P)^(-P)>* = ^ "^-^ • ( 3 - 19 ) 

This produces a Stingl propagator which has been observed before, [18J, and which will be a 
common feature in other sets of propagators. However, the gluon propagator is suppressed 
similarly to the original Gribov propagator. For the Bose ghost there are several massless poles. 
It is these features of suppression and location of massless poles which is the central focus of 
this propagator analysis. If instead we include the mass terms already proposed in [36l E7J EH] 
we reproduce those results but record them in our current notation for completeness. We have 



(A°(p)Al(- 


-p))q = 


(A a ^ui- 


-p))q = 


(A a »(p)pl c (- 


-p))q = 


<C 6 (P)C(- 


-p))q = 


{^ b (p)pl d (- 


-p))q = 




-p))q = 



[(P 2 ) 2 + PqP 2 + °A1 

2 jabc^2 



pAp) 





^ac^bd jabe. jcde^A 

^ + r„2 , ,.2ir^2A2 i ..2„2 i 7T~Zi'- ? V(p) 



b 2 + 4r [P 2 + P 2 Q ][(P 2 ) 2 + P 2 Q P 2 + C A ^] 



Xacsbd 

«(pK d (-p)) Q = - ^ + fl 2 ] ^u ■ (3.20) 

Here the gluon propagator does not vanish at zero momentum and there are no massless poles 
in any of the colour channels. It was in part this gluon propagator freezing which justified 
examining the condensation of the associated BRST invariant operator originally, |36[ [37[ EH] ■ 
As the expressions for both the separate W and S channel propagators are cumbersome we have 
recorded them in Appendix B for the case of SU(N C ). However, the case of the 1Z channel is 
similar to that for Q since 



[(P 2 ) 2 + P 2 n p 2 + C A1 4 

if abc -/ 2 

[(P 2 ) 2 +74F+cu7 4 r 

(A^p)p^(-p)) n = 

(tabbed, « ^S bd f ahe f cde [p 2 n P 2 + C A1 '\ 



(Al(p)£ c {-P))n = MU V' X n^ p v><?) 
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f abe f cde 2 

+ c AP 2 [ P 2 + p 2 n ] L ^ [p) 

(Cf(p)f^(rP))n = 

sacxbd fabe fcde..2 

(pf(p)pl d (-p))n = K b ( P )rt d (-p))K = ~ —s-ri^ + r 2/2, 21 V (3-21) 

P ^AP [P 1 p-ji\ 

and we note that the inclusion of the factor of C A in (|3.ip eases comparison. Like (|3.20p the 
gluon propagator freezes. However, the main difference is that there are massless poles in, for 
instance, the propagator. Whilst it was these massless poles which became enhanced when 
the gap equation was satisfied by 7 in the pure Gribov-Zwanziger Lagrangian we will show later 
that there is no similar enhancement in this case. The situation for the V channel has parallels 
with the previous set since 



raft 2 

{A-(p)At(-p))v = - [(p 2 ) 2 + C A ^ {p) 



{Al{ V )pl c {-p)) v = 

Sjacxbd f abe t cde ^4 sab zed. .2 

(tt(p)tf(-p))v = o 

Xacxbd xab£cd..2 

(Pl b (p)pt d (-P))v = K b (p)ut d (-p))v = - + NAp 2 [v 2 + ^ (3-22) 

producing massless poles but gluon suppression instead of freezing. By contrast the T channel 
has gluon freezing but no massless poles because 



(A»( P )Al(- 


-p))t = 


{Al{p)&{- 


-p))t = 


(A a »(p)pl c (- 


-p))t = 




-p))t = 



b ah \p 2 -p 2 T \ 



pAp) 



[(p 2 ) 2 -^P 2 + CA7 4 ] "" 

j fabc^2 



[(p2)2_ fl 2 Tp 2 + CAl A ] 





pAp) 



^ac^bdp2 jabe jcde^i 

i 7 ?^ + O ..2 ir/„2A2 ,.2 „2 i FTTZT-'VCp) 



W) 2 -p\\^ \p 2 -p 2 t ][(p 2 ) 2 -p 2 t p 2 + c a ^\ 



b ad b bc p l T 

(tf(pK d (-p)h = o 

ftaczbd zadzbc 2 

(p a » b (pK d (-p)h = (^{p)ul d {-p)) T = ~ + [(p2)2 _ J 4] ^ • (3-23) 

Though, of the sets we have recorded this appears to be the one which is least likely to be 
realistic since a type of tachyonic pole appears which would violate causality. 

As there is interest in effective gluon masses we record the propagators for two more general 
situations. Those with a non-zero p? x and either a non-zero p 2 Q or a non-zero p 2 ^. The former 
case has been discussed in [36] but as we will concentrate on the 1Z channel we will examine the 
consequences of a gluon mass for this case briefly as well. The motivation for a non-zero p 2 x 
originates from the same argument as that leading to (I2.16P . If we integrate the gluon propagator 



12 



then one obtains a non-zero vacuum expectation value for ^A^A " 11 . For completeness and to 
allow us to contrast with other situations we note that the XQ case is 

§a,btp2 _|_ ^2 j 



(K(p) P b u c (-p)) XQ = 

[p 2 + p 2 q ] ' 



ftacsbd 

(C b (p)C d (-p))xQ = - ,J , „ 21 ^ 



fabe fcde^i 

\P 2 + P 2 Q ] [(P 2 ) 2 + + ^ 2 S )P 2 + p\p\ + c Al ^\ v W 



<£?(p)^(-p))*Q = 

xacxbd 

(P; b (p)p c u d (-P))XQ = (rfW?{-p))xQ = ~ [p 2 + fl 2 Q] V^ ■ (3-24) 

The presence of the non-zero p x does not alter the properties significantly from the pure Q 
channel case. There is still gluon freezing and no massless poles in the full set. The structure 
of the Stingl propagator is unsurprisingly affected. Equally when p 2 x and p 2 ^ are both non-zero 
the propagators behave in essence in the same way as the 7Z channel ones since 



ir*M?*(- n \\ 5<1CSbd n f abe f cde [p 2 nP 2 + »Wn + C A ^] 



fabe fade 2 
, J I P"R. j (r \ 

(i ab (p)p cd (-p))xii = 

racrbd sabefcde.,2 
(P; b (p)p c u d (-P))XK = {ut{p)Q?{-p))xK = - —V»» + C ^ 2+ ^] VM • 

(3.25) 

Clearly there is gluon suppression together with massless poles. For both (|3.24p and f|3.25|) one 
could have the situation where if p 2 x = — p 2 Q or p x = — p 2 ^ then one would return to a Gribov 
type denominator in the gluon propagator. However, the main point to appreciate from this 
analysis is that there is another possibility of obtaining a frozen gluon via the condensation of a 
BRST invariant dimension two operator which has a different colour structure to that considered 
in [361 1371 [38] . This derives from (|3.2ip and we will refer to it as the 1Z channel mechanism for 
the moment and study it in more detail in the next section. 

Prior to that we need to briefly discuss the renormalization of the general operator Q abcd 
since we will be performing loop computations. It has already been shown that the contraction 
with 5 ac 5 bd produces a renormalizable operator, [361 1371 [38] . However, the more general oper- 
ator is clearly renormalizable multiplicatively in the Landau gauge by the same reasoning. To 
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partly verify this and check our internal conventions we have renormalized O to two loops 
by inserting the operator in either a t/>° b or a io^ b 2-point function where the momentum flows 
in through one of the legs and out through the operator itself. Our momentum configuration 
allows one to apply the Mincer algorithm of [H] but recoded, [49], in the symbolic manipu- 
lation language Form, [50]. As the two external fields in the Green's function each carry one 
Lorentz index then we have to contract all the Feynman integrals with the tensor rj^/d so that 
Mincer can be applied to a Lorentz scalar. The factor of d in the denominator arises from the 
normalization. Concerning the colour indices we do no projection on the four free indices of 
the operators which allows us to check that the divergent part of the Green's function is only 
associated with the colour structure of the original Feynman rule for the operator. In essence 
this is how we check that the renormalization of the general operator is multiplicative. For the 
Green's function with § ab legs there are 4 one loop and 160 two loop Feynman diagrams. The 
respective figures for the Green's function are 1 and 25. The diagrams are generated by the 
Qgraf package, [51] . before being converted into Form input notation. Then the Form ver- 
sion of the Mincer algorithm, |49j . is applied. As we are using an automatic Feynman diagram 
computation procedure we use the method of rescaling of bare quantities such as the coupling 
constant in order to deduce the overall final renormalization constant for the Green's functions 
we renormalize, [52] . For both Green's functions we find that the operator is multiplicatively 
renormalized and with the same renormalization constant Zq. Not unexpectedly like the earlier 
colour contraction the more general operator Q abcd satisfies the Slavnov- Taylor identity 

Z^Zo = 1 (3.26) 

to two loops in MS where Z^ is the renormalization constant of the original localizing ghost 
field § ab . The renormalization constants for the real and imaginary parts of this field § ab are 
equivalent. 



4 1Z channel mass. 



Having considered the different forms the propagators of the spin-1 sector can take when there 
are a variety of single mass terms originating from the BRST invariant operator Q abcd ^ we now 
focus on the 1Z sector mass in detail. Like the Q case it has a frozen gluon propagator. The 
aim is to give evidence that there is a dynamical origin for such a mass as the previous analysis 
merely assumed the existence of such additional mass terms. The approach we take is the same 
as that of [37] for the Q channel. Indeed given the similarities between the Q and 1Z colour 
contractions of O abcd a large amount of the results of [37] can immediately be transferred to 
the present case without detailed re-analysis. For instance, the inclusion of the general operator 
(|3.ip does not destroy the renormalizability of the pure Gribov-Zwanziger Lagrangian, (|2.4p . 
Moreover, we use the same procedure to examine the dynamical origin of an 1Z channel mass 
which is the construction of an effective potential for the operator. From this it will transpire 
that there is a non-zero vacuum expectation value for the operator which therefore condenses 
to produce the mass term we analysed previously. However, it turns out that from the way we 
have set up the 1Z term of (|3.2p that the effective potential is formally the same as that for the 
Q case whence we can merely translate all those results to the present case. For instance, the 
starting point of the application of the local composite operator formalism, [55] HU] HI] [22], to 
the operator is the one loop energy functional W[J] where J is the source coupling to the 7Z 
channel operator. Summing up the conventional set of contributing one loop diagrams produces 



W[J] 



dN Al 4 

2g 2 



+ 



dN A (Jj 2 



+ 



(d - 1)N A 



[ In 


k 2 


Ik 





k 2 + 



Ca1 a 
[k 2 + J] 



+ 0(g 2 ) (4.1) 
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where £ is the local composite operator parameter, [39ll4C H l4"T | l42]. Clearly this is formally similar 
to that for the Q channel and hence we merely recall the subsequent properties. Defining 

•« - m 

then the effective action, T[a], is constructed from the Legendre transformation 

T[a] = W[J] - J d i xJ(x)a(x) . (4.3) 

Setting 

a{x) = cr + a(x) (4.4) 

where 

(To = 5 ( 4 - 5 ) 

then we find the same non-zero value for the condensate as before, 1371 . 



a 



64vr 



(4.6) 



With this observation then the corresponding subsequent analysis of [37] in respect of the gluon 
propagator will hold for the 1Z channel mass. This is due to the fact that the gluon propagators 
of (|3.20p and (|3.2ip are formally the same. Therefore, the estimate for the value where the gluon 
propagator freezes, which is in qualitative agreement with the lattice, applies in this case too. 

Given this complete parallel between the two cases it might appear that there is no justifica- 
tion in posing an alternative way of having a frozen gluon in a refinement of the Gribov-Zwanziger 
formalism. However, there is a key difference and it resides in the Bose ghost sector. From (I3.20P 
and (|3.2ip there are massless poles in the £^ propagator for the latter case but not for the former. 
Recently, Zwanziger has argued, |43| . that in the pure case there is an enhancement of the Bose 
ghost in the infrared which is a non-perturbative property of the theory. The argument is based 
on the spontaneous breaking of the BRST symmetry in a theory where fields are constrained 
by the horizon condition. The structure of the propagator has been examined at one loop in 
the MS scheme in the zero momentum limit and an enhanced emerges, [21]. There is also 



enhancement in the case too. In essence the starting point for this is the massless poles of the 
original propagators, (|2.14p . However, one needs to compute the one loop corrections to all the 
2-point functions of the fields in the zero momentum limit and apply the gap equation satisfied 
by 7. Then the leading momentum term vanishes to ensure that the resultant propagators have 
a dipole behaviour in the infrared as opposed to the canonical behaviour of a massless field. The 
analysis for is hampered by the colour tensor structure of the 2-point functions and hence 
only certain colour channels of its propagator enhance. 

Therefore, given that the structure of the propagators in (I3.20P and (13.2ip are different 
it is worth examining the structure of the 7Z propagators in the infrared limit. However, as the 
gap equation is important for this we need to compute the one loop expression for 7. This is 
achieved by evaluating (|2.9|) using the mixed propagator of (|3.2ip . As the denominator factors 
are of a Stingl form in order to evaluate the basic Feynman integrals we need to first rewrite 
the factor formally in a more conventional fashion. Therefore, we define two additional mass 
parameters, given by 

/4 = 4 + M- , C A1 A = 4/i 2 _ . (4.7) 
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This choice is motivated by the relation 



Solving the relations produces the mapping 



(4.f 



/4 



/4 



/4 " ± C A1 A 



(4.9) 



With this factorization it is a straightforward matter to determine the MS gap equation and we 
find 



C A 



5 In 



Ca7 4 
M 4 



3/4 



/4-4CU 7 4 



:hl 



/it 



lit 



a + 0(a 2 



(4.10) 



In the limit of /4 — )• one recovers the gap equation for 7 of the pure Gribov-Zwanziger case. 
This expression is formally the same as that for the Q channel, [36} 137] . if one simply replaces 
/4 by ii 2 q. This is not unexpected since the mixed propagators are formally the same at this 
order. At two loops one would expect the gap equations to be different merely because at that 
order the other spin-1 propagators will appear in the Feynman diagrams. 

Before considering the infrared structure of the propagator we first examine the Faddeev- 
Popov ghost propagator. This is partly to highlight the technique one follows but in a case which 
is not complicated by the colour tensor structure as well as to verify the loss of ghost enhancement 
which has to be a feature of the 1Z channel in order to be consistent with the evidence from the 
lattice. However, the situation for the behaviour of the Faddeev-Popov ghost propagator in the 
infrared is the same as that for Q since at one loop only one diagram contributes to the 2-point 
function. Defining the ghost propagator in terms of its form factor, D c (p 2 ), by 



Dc(p 2 

p2 



xtib 



(4.11) 



we have 



D c (p 2 ) 



1 - C A 



+ 



In 
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Ca1 A 

(P 2 ) 2 
Ca1 a 

1 



+ 



3/4 



:lll 



lit 



CaY 



In 



4^7/4 -4C A7 4 



p 2 + 0((p 



„2\2 



+ 0(a 2 ) 



(4.12) 



in the limit as p 2 — > where we have included the 0(p 2 ) contribution. Ordinarily this term 
would govern the infrared behaviour of the propagator in the infrared as the leading order term 
in momentum would vanish when 7 satisfies the gap equation. That does not happen due to 
the O(iift) term which has the opposite sign to the analogous term which appears in the gap 
equation, (|4.10p . Therefore, as with the Q channel there is no Faddeev-Popov ghost enhancement 
for 1Z either. Again this is in keeping with the lattice data which sees a minor variation from 1/p 2 
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behaviour of the ghost propagator in the infrared, [231 [Ml E3 E3 EZl ESI ESI [30], Ell As the 

colour structure of the p ab and io^ b propagators is equally as trivial as that of the Faddeev- Popov 
ghost then the propagators of these fields are not enhanced either. 

It is interesting to trace the origin of the sign discrepancy that prevents the enhancement 
happening at one loop. For instance, the ghost 2-point function correction is purely from one 
Feynman diagram and that includes a gluon propagator. By contrast evaluating the horizon 
vacuum expectation value involves the mixed propagator at one loop. Therefore, from examining 
these propagators it transpires that it is the part of the gluon propagator leading to its freezing 
which is responsible for the difference in the signs of the respective terms in the gap equation and 
the ghost 2-point function. In other words if the /4 term in the gluon propagator numerator was 
absent then there would be one loop ghost enhancement even with the Stingl type propagator 
whose denominator has to be the same in both the gluon and mixed propagator. With this 
observation it is simple to determine which of the various channels we have introduced in an 
earlier section will contain enhanced Faddeev-Popov propagators. This will happen for X, V 
and S. So, interestingly, a massive gluon but with a Gribov width would be infrared suppressed 
whilst satisfying the Kugo-Ojima confinement criterion at one loop in this context. Whether 
there is enhancement in various colour channels for the fields and if so which ones, is more 
complicated to determine since it would require the structure of the 2-point function as well 
as all the other 2-point functions. 

To illustrate the complexity of such an exercise even in a simple case we return to the 1Z 
channel and record the explicit one loop form factors for all the spin-1 field 2-point functions. 
We have 
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where /2 is the mass scale introduced to ensure that the coupling constant remains dimensionless 
in d-dimensions as we are using dimensional regularization. The divergences have been removed 
by the MS scheme prescription. As expected is similar to the Faddeev-Popov ghost 2- 
point function. In the pure case this function was solely responsible for producing the overall 
enhancement of the propagator observed in |21| which was consistent with the general analysis 
of [13]. However, since the leading term of the momentum expansion does not vanish when the 
gap equation is satisfied this leads to the absence of enhancement for in the TZ channel. To 
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be more specific we find the leading behaviour in the p 2 
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which is derived using the standard procedure from the coefficient of the leading term in the zero 
momentum limit after the gap equation has been set. It is the same as that for the Faddeev- 
Popov ghost. As expected there is no enhancement and unlike the pure Gribov-Zwanziger 
case the leading momentum behaviour involves a logarithm of the momentum. However, the 
Bose ghost propagators diverge in the same way as the Faddeev-Popov ghost. Interestingly the 
colour tensor structure is the same as that for the enhanced £^ propagator of the pure Gribov- 
Zwanziger case discussed in 
that case is 



[2T] and we recall the leading behaviour of both Bose ghosts in 
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So far in our discussions we have noted that the properties of the Q and 1Z channel masses are 
the same. However, (|4.14p clearly represents the first departure from that similarity. For the Q 
case all colour channels of the propagator freeze to finite values in the zero momentum limit. 
By contrast the 1Z channel propagator neither freezes nor enhances. Instead its behaviour is 
in essence the same as that of the Faddeev-Popov ghost since the coefficient of the leading term 
is the same but differs in the colour tensor structure. That the same tensor structure as (14.161) 
emerges at leading order for both cases is a consequence of the peculiarities of the inversion of 
the matrix of colour structures. However, the adjoint colour projection that Zwanziger focused 
on in [33] will not enhance nor diverge but freeze to an non-zero value. 

There are several main consequences of our 1Z channel analysis. First, it appears that there 
are now two possibilities of modelling lattice data by refining the Gribov-Zwanziger Lagrangian 
with one additional mass operator. Both the Q and 1Z channels reproduce gluon freezing to a 
non-zero value and ghost non-enhancement. However, to determine which one is correct would 
require a lattice computation of the ^ propagator in the zero momentum limit. It would seem 
to us that this would be a non-trivial exercise since the lattice gauge fixing procedure used 
pays no attention to a field, which only appears in the Gribov-Zwanziger Lagrangian and 
is necessary to localize the non-local horizon operator. However, one could consider instead the 
correlation of a non-local projection of the gluon field itself since the equation of motion from 
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produces the relation 
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(4.18) 



where we have been careful in including the colour indices of the inverse Faddeev-Popov operator 
to correct an error in earlier work, [20, 21\. Again this would appear to open up other difficulties 



20 



since the presence of a non-locality via the Faddeev-Popov operator may be hard to define in 
a discrete spacetime in an unambiguous way. Aside from that it may even be restrictive both 
financially and computationally to produce accurate enough data in order to determine the 
behaviour at zero momentum definitively. Currently, the state of the art to obtain the gluon 
propagator behaviour at zero momentum involves a new formulation of the lattice definition 
of the linear covariant gauge, [53]. Another problem relates to whether (|4.18p is the proper 
definition of the field whose correlator is the relevant object to examine. For instance, (14.181) 
assumes there is no mass term initially for the field and that the mass term has a dynamical 
origin. However, if one has the mass term present then it could be argued that (|4.17p should be 
replaced by something such as 

(d v D u ^) ab + ^r fee / c *^ = rffahc^ (4 19) 

Though in the infrared limit this would effectively be the same as an adjoint projection of the 
gluon. So it is not clear in this case whether this would allow one to make a clear statement on 
the behaviour of the propagators we have considered, {^{p)H d {—p))i and {p a ^(p)p cd {—p))i- 



5 Discussion. 



There are several main features arising out of our analysis which was to attempt to model 
lattice data with a more natural perturbation of the original Gribov-Zwanziger Lagrangian with 
a different BRST invariant dimension two operator compared to |36l [37], [38] • First, if one is to 
reconcile the lattice and Schwinger-Dyson observations that the Landau gauge gluon propagator 
freezes to a non-zero value and the Faddeev-Popov ghost does not enhance, by modelling with 
a refined Gribov-Zwanziger Lagrangian, then it turns out that there is more than one way to 
do this. The original analysis of [361 EH EH] used only one colour projection but it has been 
shown here that the alternative 1Z channel projection can equally accommodate a non-zero frozen 
gluon propagator and unenhanced ghost propagator. If one wished to determine which of these 
single operator extensions was consistent with the lattice results then the test resides in the 
infrared behaviour of the propagators of the Bose localizing ghost. In the earlier Q channel 
solution the £J 6 propagator freezes to a non-zero value. By contrast in the 7Z channel case the 
propagator of this field has the same momentum behaviour in the infrared as the Faddeev-Popov 
ghost propagator. Given this it might be better in future to refer instead to the Q channel as 
the gluon-like massive or decoupling scenario and that of 1Z as the Faddeev-Popov ghost-like 
massive or decoupling solution. In some sense the resolution by numerical work may not in fact 
be possible given the amount of computing resources which would be required to determine the 
correlation of a non-local projection of the gluon field. Moreover, both sets of localizing ghosts 
are an inherent feature of accommodating the original non-local Gribov horizon operator and 
as such would have no parallel or direct concept in a lattice construction. 

We have also noted that the inclusion of the dimension two gluon mass operator, which is 
BRST invariant as well, into the propagators for each of the Q and 1Z channels does not alter 
their behaviour from the situation when the gluon operator was absent. This in fact opened 
a wider question as to whether all the possible colour projections should not be considered 
simultaneously. As the gluon and Faddeev-Popov ghost propagators are the only two which 
have been analysed numerically the fact that the gluon freezes to a non-zero value can also be 
accommodated by contributions from four different colour channels as can be seen in Appendix 
A. From a numerical point of view the freezing to a non-zero value for the gluon propagator 
can never resolve each of the four different mass scales which arise there. Again only data 
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on the remaining spin-1 propagators could ever possibly determine this. It would require, for 
instance, a substantial amount of numerical fitting of all the mass parameters. Though as 
there are more than seven different propagator form factors then this over-redundancy would 
provide independent consistency checks on the seven mass parameters where we include the 
gluon mass parameter in the counting here. However, we should recall one of our underlying 
assumptions in this context. In (|3.ip we have taken the masses for p° 6 , and u^ b to be the 
same in order to ensure our additional operator is BRST invariant. If the BRST symmetry 
is broken then there is no reason why, for example, any of the masses of the localizing fields 
should be equal. If this were the case then the propagators we have discussed throughout 
would be much more complicated. Moreover, this would become apparent in all the different 
propagator form factors analogous to those of Appendix A. Again to test this scenario would 
appear to be computationally impractical for reasons we have already mentioned. However, one 
alternative theoretical way to shed some light on the interplay of the different colour projections 
would be to extend the effective potential analysis of the single projection case for the Q and 
1Z channels to include all seven cases simultaneously. Whilst it is a non-trivial task to apply 
the local composite operator method for this, finding a stable absolute minimum of the effective 
potential would indicate which is the most energetically favoured solution or solutions. This is 
currently in progress. However, it seems from earlier experience, with the construction of the 
enhanced Bose ghost propagator in the pure Gribov-Zwanziger Lagrangian, that the subtleties 
of the intricate nature of the general colour tensors of the localizing ghosts have a significant 
effect on the infrared structure of this Lagrangian. One only has to recall the decompositions 
(|3,15p to see how the effective potential construction will mix up colour tensors of the operator 

Qabcd 

Finally, in providing a complete analysis for each of the seven individual colour projections 
separately we have noticed several features in the interplay of a frozen gluon propagator and a 
non-enhanced Faddeev-Popov ghost propagator with the structure of the original propagators 
when the appropriate gap equation is implemented at one loop. Specifically if the propagators 
of the spin-1 fields have a massless pole in a colour or Lorentz channel then whether it becomes 
enhanced in the infrared depends on whether the gluon propagator freezes to a non-zero value 
or not. In the case of the original Gribov-Zwanziger Lagrangian this was the case as is now 
evident in the recent work of |43} I21j for the localizing Bose ghost. Here this particular feature 
emerges in several of the cases. For instance, with a gluon mass operator only or in the S 
channel only then when the appropriate horizon condition for 7 is satisfied there appears to 
be ghost enhancement. Whilst there is currently no lattice evidence for this scaling scenario 
since the massive or decoupling solution appears to be favoured by many different simulations, 
^|21l|23|2g|2Zl|2Sl|2g|^EBI321ISll531, itis perhaps worth noting that if instead a scaling 
solution had been found then for this situation the pure Gribov-Zwanziger structure may not 
actually have been the unique explanation. 



A Full propagators for 577(3). 

In this appendix we record the explicit expressions for the propagators for the specific colour 
group SU(3) for all channels together. First, the transverse sector is 
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Clearly there are no massless poles provided the combinations of the various masses in each of 
the denominators do not accidentally sum to zero. Further, the gluon propagator is transverse 
and freezes to a non-zero value in the infrared when again there are no cancellations between 
the parameters. 



B W and S channel propagators for SU(N C ). 



As the explicit expressions for the propagators for each of the W and S channel cases are 
complicated for an arbitrary colour group we present the expressions for SU (N c ) only here. 
First, the W propagators are 
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Clearly the gluon propagator freezes to a non-zero value and the t\ ab and p ab propagators have 
massless poles in various colour channels. More specifically the dominant part of each propagator 
as p 2 — > is 
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Interestingly the colour tensor structure of the dominant infrared behaviour of each propagator 
is the same as that for the enhanced propagator in the pure Gribov-Zwanziger theory after 
the gap equation has been implemented. Though the Lorentz structure differs. As the gluon 
propagator is essentially the same as that in the 1Z channel case then the same features will arise 
when the one loop gap equation is computed for this case. In other words the massless poles 
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in the full set of W channel propagators will not enhance for similar reasons to those of the 1Z 
case. Finally, just to be complete for this channel we note that the leading order behaviour of 
the adjoint colour projected ^ propagator as p 2 — > is 
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Clearly the transverse part of this particular correlator vanishes in the infrared but the longitu- 
dinal part freezes. 

For the S channel the propagators are equally as involved since 
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By contrast there is no gluon freezing but the massless poles are distributed to the same colour 
channels as W. However, the colour tensor structure for the propagator at leading order in 
the zero momentum limit is more akin to that observed for the enhanced propagator in the 
pure Gribov-Zwanziger case. Specifically as p 2 — > we find that 
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where we include that for p ab for completeness. Whilst the same colour tensor structure emerges, 
the associated enhancement can only be determined when the one loop corrections are computed 
to all the 2-point functions. If we examine the colour adjoint projection of the propagator 
here then we find that at leading order in the zero momentum limit 



(f abp C b (p)f cdq C d (-p))< 
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So the transverse part is finite and indeed vanishes as p 2 — > whilst a massless pole is present 
in the longitudinal part. As discussed earlier we would expect that this will be enhanced when 
the loop corrections are included in the 2-point functions since here the gluon propagator is 
suppressed. 
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